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We report measurements of B-meson decays into two- and three-body final states containing two 
charmed baryons using a sample of 230 million T(4S) — > BB decays. We find significant signals in 
two modes, measuring branching fractions B(B~ -> At A~ K~) = (1.14 ± 0.15 ± 0.17 ± 0.60) x 10~ 3 
and B(B~ -» S°AJ) x B(S° -> S~tt + ) = (2.08±0.65±0.29±0.54) x 10" 5 , where the uncertainties 
are statistical, systematic, and from the branching fraction B(A~£ — > pK~n + ), respectively. We 
also set upper limits at the 90% confidence level on two other modes: B(B° — > S£ A^) x B(E^ — > 
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5 7T+7T+) < 5.6 x 10 5 and B(B° -> A+ A~ K°) < 1.5 x 1CT 3 . We observe structure centered at an 
invariant mass of 2.93 GeV/c 2 in the A^K~ mass distribution of the decay B~ — > At A^ K~ . 

PACS numbers: 13.25.Hw, 12.15.Hh, 11.30.Er 



Bottom (B) mesons are heavy enough to decay into 
charmed baryons, and do so at a rate of roughly 5% [HEl. 
The dominant decay mechanism is via b — > W~ transi- 
tions, with W~ coupling to cs oiud 0], both of which are 
Cabibbo-allowed. Theoretical predictions for the branch- 
ing fractions of B mesons to baryon-antibaryon pairs 
have been made within the diquark model Q and with 
QCD sum rules These suggest that decays to two 
charmed baryons (B — > X C \X C 2) and to one charmed 
baryon and one light baryon (B — > X C \X2) have compa- 
rable branching fractions, of the order of 1CP 3 for indi- 
vidual modes. 

Several inclusive measurements of i?-meson decays to 
charmed baryons have been made In particular, 

the BABAR Collaboration recently performed an inclusive 
analysis of A^ production in which flavor tag information 
was used to identify whether the came from a B or 
a B meson @ . It was found that about a third of all A+ 
were from B mesons with anti-correlated flavor content 
(i.e. b — > c rather than b — > c transitions), consistent 
with a substantial rate of b — > gs decays. Inclusive stud- 
ies of the and A£ momentum spectrum [3, H @| also 
support a substantial rate of baryonic b — > ccs decays 
such as B~ — > 5°A~. However, inclusive studies alone 
cannot fully establish this, since the momentum distribu- 
tions can also be reproduced with carefully tuned sums 
of b — ► cud processes. Therefore, exclusive measurements 
are needed. These require very large samples of i?-meson 
decays and have only recently become feasible. 

The Belle Collaboration has reported results on B de- 
cays to final states with two charmed baryons in both 
two- and three-body modes [H, They measured 
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± 1.1 ± 



B(B~ -> 5°A~) x_B(~° -> H-7T+) = (4.o_ . 9 
1.2) x 10- 5 and B(B° -> S+Az) x B(5+ -► ~-ir+ir+) = 
(9.3± 3 2 l ± 1.9 ± 2.4) x 1CT 5 Q. Assuming that B(S° -> 
S~tt + ) and B(3+ — > E~n + n + ) are of the order of 1%- 
2% [11[ , these results are compatible with the prediction 
that B(B~ -> S°A~) and B(B° -> S+Az) are ©(lO" 3 ). 
This is in stark contrast to the branching fractions of 
singly charmed decays, such as that of B° — ► A+p which 
is (2.2 ± 0.8) x 10 -5 , smaller by two orders of magni- 
tude [12]. The branching fractions of the three-body 
processes B — ► A+A~K were also found to be large: 
B(B~ -> A+A-K-) = (0.65ig;Jg ± 0.11 ± 0.34) x 10" 3 
£(B° yl+i-i?°) = (0.79±g;||±0.12±0.42) x 10" 3 ftjfr. 
Explanations for these widely varying values have been 
proposed [H, Q ■ It was suggested that a kinematic sup- 
pression may apply to decays in which the two baryons 
have high relative momentum, since this requires the 
exchange of two high- momentum gluons. The rate of 



B — > A+A~K decays could also be enhanced by final- 
state interactions, or by intermediate charmonium reso- 
nances. 

In this paper, we present measurements of the branch- 
ing fraction of the decays B~ — > A+A~K~, B~ — > 
E°A-, B° -> E+A~, and B° -> A+AzK , and inves- 
tigate three-body decays for the possible presence of in- 
termediate resonances. The data were collected with the 
BABAR detector [l5[ at the PEP-II asymmetric-energy 
e + e~ storage rings and represent an integrated lumi- 
nosity of approximately 210 fb" 1 collected at a ccnter- 
of-mass energy y/s = 10.58 GeV, corresponding to the 
mass of the T(4S) resonance. The BABAR detector is 
a magnetic spectrometer with 92% solid angle tracking 
coverage in the center-of-mass frame. Charged parti- 
cles are detected and their momenta are measured in a 
five-layer double-sided silicon vertex tracker and a forty- 
layer drift chamber, both operating in a 1.5 T mag- 
netic field. Charged particle identification (PID) is pro- 
vided by the average energy loss (dE/dx) in the track- 
ing devices and by an internally reflecting ring-imaging 
Cherenkov detector. Photons are detected with a CsI(Tl) 
electromagnetic calorimeter. The instrumented flux re- 
turn for the solcnoidal magnet provides muon identifi- 
cation. Simulated events with B mesons decaying into 
the relevant final states are generated with EvtGen [l6| 
and PYTHIA (T3] , while GEANT4 is used to simulate 
the detector response. Inclusive Monte Carlo (MC) sam- 
ples of T(45*) and e + e~ — ► qq (q = u, d, s, c) events at 
^/s = 10.58 GeV are also used, corresponding to more 
than 1.5 times the integrated luminosity of the data. 

The candidates are reconstructed in the three de- 
cay modes pK~ir + , pK®, and Air + ; candidates in the 
two decay modes S~tt + and AK~ir + ; and E+ candidates 
in the decay mode Ez~it + tt + . We begin by reconstruct- 
ing the long-lived strange hadrons: Kg — > 7r + 7r~ and 
A — > pir~ candidates are reconstructed from two oppo- 
sitely charged tracks, and S~ — > Air~ from a A candidate 
and a negatively charged track. In each case, we fit the 
daughters to a common vertex and compute their invari- 
ant mass. The mass is required to be within 3er of the 
central value, where a is the experimental resolution and 
is approximately 4.0, 4.5, and 6.0 MeV/c 2 for A, and 
E~ , respectively. Candidates with a x 2 probability be- 
low 10 -4 are rejected. For A candidates, we also require 
the daughter proton to satisfy PID criteria. The mass of 
the K®A, or E~ candidate is constrained to its nominal 
value [l| for subsequent fits. 

We suppress background by requiring the transverse 
displacement between the event and decay vertices to be 
greater than 0.2 centimeters for K®, A, and S~, each 



5 



V/c 


35 


-> 




s 


30 






— 


25 






s 


20 


> 




H 


15 




10 



0.2 
0.15 

0.1 
0.05 
-0 
0.05 
-0.1 
0.15 
-0.2 c 



,(a) 

Ut„i i i iiUkiji \ 






i 


~~ (b) 


it i 



wc compute the energy-substituted mass toes 

2 



(a/4- 

p^") 1 / 2 and the energy difference AE = E* B - y/s/2, 



'5.2 5.22 



5.24 5.26 5.28 
m ES (GeV/c 2 ) 



20 40 60 

Events / (8 MeV) 



FIG. 1: The m E s and AE distributions for B~ -> A+ A~ K~ 
candidates, summing over five different final states. Plot (b) 
shows the scatterplot of tties vs. AE, and (a) and (c) show 
the mES and AE projections for \AE\ < 0.022 GeV and for 
mES > 5.27 GeV/c 2 , respectively. The dashed horizontal and 
vertical lines in (b) indicate the signal regions used for the 
projection in (a) and (c), respectively. 



of which travels several centimeters on average. We also 
require that the scalar product of the displacement and 
momentum vectors of each hadron be greater than zero, 
and that the transverse component of the displacement 
vector of a S~ candidate be smaller than that of its A 
daughter. 

Next, we reconstruct the charmed baryons A+, S" n , and 
in the decay modes listed previously. In each case, 
we fit their daughters to a common vertex, require the 
invariant mass of the charmed baryon candidate to be 
within 18MeV/c 2 (approximately three times the exper- 
imental resolution) of the nominal mass [![, reject can- 
didates with a x 2 probability below 10~ 4 , and then con- 
strain the masses to their nominal values. Wc also require 
that daughter kaons and protons of the charmed baryons 
satisfy the PID criteria for that hypothesis. 

We reconstruct B-meson candidates in the following 
final states: A+A~K', S°A~ , S+A~, and A+A~K°, 
fitting the daughters to a common vertex and requir- 
ing that the x 2 probability is at least 1CP 4 . We also 
apply the kinematic and PID requirements mentioned 
above to the K° and K~ daughters of the B mesons. 
Because the branching fraction and efficiency are higher 
for A+ — > pK~n + than for the other A+ decay modes, we 
use only final states in which at least one A£ or A~ de- 
cays to pK~n + or pK + n~ . For each B-meson candidate, 



where p* B , E B , and y/s are the momentum and energy 
of the B meson and the e + e~ collision energy, respec- 
tively, all calculated in the e + e~ center-of-mass frame. 
For a correctly reconstructed signal decay, the toes dis- 
tribution peaks near the nominal mass of the B meson 
with a resolution of approximately 2.5MeV/c 2 , and AE 
peaks near zero with a resolution of 6.0-7.8 MeV depend- 
ing on the final state. Figure [T] shows the Toes and AE 
distributions for B~ — > A+A~K~ candidates. 

Background arises from several sources, including mis- 
reconstructed B decays to two charmed baryons, B de- 
cays to a single charmed baryon, e + e~ — > cc events con- 
taining charmed baryons, and random combinations of 
tracks. We use inclusive MC simulations and events from 
the sidebands of toes, AE, and charmed baryon mass 
in data to study the background. We consider as back- 
ground i?-meson decays with the same final state that 
do not proceed via an intermediate charmed baryon — 
for example, B~ — > E®pK + ir~ misinterpreted as B~ — > 
E®A~. Decays of this kind are distributed as signal in 
Toes and AE but have a smooth distribution for the mass 
spectrum of the misreconstructed charmed baryon, un- 
like signal decays which also peak in the charmed baryon 
mass. In studies of the S c and A c mass sidebands, we 
find no evidence for these processes and conclude that 
their contribution is negligible. 

Another important source of background is feed-down 
from related processes. The B meson can undergo a 
quasi-two-body decay via an excited charmed baryon 
such as B — > S*A~, or a non-resonant multi-body decay 
such as B — > E c A~-k. These events have similar distri- 
butions to the signal for toes and the charmed baryon 
invariant masses, but are displaced in AE by an amount 
that depends on the final state but is generally more than 
50 MeV. We remove these backgrounds by requiring that 
signal candidates satisfy lAi^l < 22 MeV. Finally, we re- 
quire 5.2 < toes < 5-3 GeV/c 2 . The average number 
of reconstructed B candidates per selected event varies 
between 1.00 and 1.14 depending on the final state. In 
events with more than one candidate, the one with the 
smallest \AE\ is chosen. We verify with MC and events 
from data sidebands that this does not introduce any bias 
in the signal extraction. Studies of simulated events show 
that l%-3% of signal events are incorrectly reconstructed 
with one or more tracks originating from the other B in 
the event; this effect is taken into account implicitly by 
the efficiency correction described later. 

The signal yields are extracted from an unbinned ex- 
tended maximum likelihood fit to the toes distribution. 
We use separate probability density functions (PDFs) for 
signal and background events. The likelihood function 
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TABLE I: Fitted signal yield, detection efficiency e, significance S, measured branching fraction B, and (for S < 2) the upper 
limit on B for each decay mode. The uncertainties on B are statistical, systematic, and the uncertainty from the branching 
fraction B(At — > pK~it + ). For final states containing S° or S+, B includes a factor of B(S° — ► S~w + ) or B(5£ — > 2~-k + -k + ), 
respectively. 



Decay Mode Signal Yield e(%) g B Upper Limit on B 



B~ -> 4+4" K~ 




74.6 ± 9.8 




9.6 


(1.14 ± 0.15 ± 0.17 ± 0.60) x 10" 


-3 


At ^ P K-n + , A~ 




42.7 ± 7.7 


7.1 


7.1 


(1.07 ± 0.19 ± 0.16 ± 0.56) x 10" 


-3 


A+ -^ P K-Tv+, A~ 




14.5 ±4.0 


8.8 


5.9 


(1.81 ± 0.50 ± 0.30 ± 0.94) x 10" 


-3 


A+^pK° s , A- -» 




11.4 ±3.7 


8.8 


4.8 


(1.42 ± 0.45 ± 0.24 ± 0.74) x 10" 


-3 


yl+ yl" 


— > 7Lt~ 


2.5 ± 1.8 


0.3 


2.0 


(0.55 ± 0.40 ± 0.09 ± 0.28) x 10" 


-3 


A+ -> Att + , A- 


P-f£T + 7r _ 


3.5 ±2.0 


6.4 


2.7 


(0.74 ± 0.43 ± 0.12 ± 0.38) x 10" 


-3 






14.0 ±4.4 




6.4 


(2.08 ± 0.65 ± 0.29 ± 0.54) x 10" 


-5 




8.0 ±2.8 


4.3 


6.1 


(2.51 ± 0.89 ± 0.29 ± 0.65) x 10" 


-5 




— > pK + ir 


6.0 ±3.4 


4.5 


2.1 


(1.70 ± 0.93 ± 0.30 ± 0.44) x 10" 


-5 






2.8 ± 2.0 


2.6 


1.8 


(1.50 ± 1.07 ± 0.20 ± 0.39) x 10" 


" 6 < 5.6 x 10" 5 @ 90% C.L. 






3.3 ±2.7 


4.4 


1.4 


(0.38 ± 0.31 ± 0.05 ± 0.20) x 10" 


3 < 1.5 x 10~ 3 @ 90% C.L. 
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FIG. 2: The fitted mES distributions observed for the de- 
cay modes (a) B~ — > At A~ , combining 5 exclusive final 
states; (b) B~ — > 2°A~, combining 2 exclusive final states; 
(c) B° -> S+A~; (d) B° -+ AtA~K°. Points with error bars 
represent the data, dashed lines the background PDF, and 
solid lines the sum of the signal and background PDFs. 



for the ./V candidates in the event sample is given by 



C 



Y[ ( n s^s("iESi) + n B VB{mvsi)) , (1) 



where S here denotes the signal and B the background, 
V is the PDF (normalized to unit integral), and n is the 
yield. The signal PDF is parameterized as a Gaussian 



function with a fixed to a value obtained from a fit to 
simulated signal events. The Gaussian mean is also fixed 
to the value obtained with simulated signal events, except 
for B~ — > At A~ K~ where there is sufficient signal in 
the data to fit this parameter. The background PDF 
is parameterized as an ARGUS function [19(. We allow 
the ARGUS shape parameter to vary within a physically 
reasonable range in the fit to the data. 

The fitted toes distributions of the four final states 
are shown in Fig. [51 Clear signals are seen in the B~ — > 
AtA~K~ and B~ — > S®A~ decay modes. A measure of 
the significance of each peak is given byS= V2Aln£ Q|, 
where A In £ is the difference in likelihood (incorporating 
the fitting systematic uncertainty) for fits where the sig- 
nal yield is allowed to vary and where it is fixed to zero, 
respectively. The results of the fits are shown in Table U 
The efficiency is determined by applying the same anal- 
ysis procedure to simulated signal events. For the three- 
body B-meson decays, the efficiency depends upon the 
distribution in the Dalitz plane. We weight the simulated 
events to reproduce the efficiency-corrected, background- 
subtracted distribution seen in data for B~ — > At A~ K~ . 
As a crosscheck, we also compute the efficiency assuming 
a phase-space distribution and find a difference of less 
than 10% in each case. 

We then obtain each branching fraction as: 



B(B^X C A-[K}) 



(2) 



where X c is the charmed baryon (At, S®, or n$j 
is the signal yield extracted from the fit to the data for 
the j th sub-mode, \\ i Bij is the product of the daugh- 
ter branching fractions, Afg is the number of neutral or 
charged B mesons, and £j is the signal detection effi- 
ciency We assume equal decay rates of the T(4S) to 
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B + B~ and B°B° Q. 

The branching fraction B(A + — ► pK~ir + ) has been 
measured previously to be (5.0 ± 1.3)% [l|. Because 
the branching fractions of and £"+ decays have not 
been determined experimentally we quote the products 
of the branching fractions, B(B~ — ► 5°A~) x B(E® — * 
E-7T+) and -> S+ylj) x B{S+ -> ~-7r+7r+). For 
the — > AK~ir + decay mode we scale the measured 
branching fraction by the ratio B(S® — ► S" _ 7r + )/S(S'° — > 
ylif-7r+) = 1.07 ±0.14 Q so that its value can also be 
expressed as the product of the same two branching frac- 
tions. 

For each decay mode, Table fl] gives the values of ns, 
e, the significance, and the branching fraction. For each 
mode with a significance below 2 standard deviations, 
we calculate the Bayesian upper limit [l[ on the branch- 
ing fraction including systematic uncertainties and ob- 
tain_6(B° -> A+A-K ) < 1.5 x 10~ 3 and B{B° -> 
S+A-) x B{E+ -> »-7r+7r + ) < 5.6 x 10~ 5 at the 90% 
confidence level. 

Table [IT] lists the main systematic uncertainties and 
their sum in quadrature. The largest uncertainty is 
from the charged track reconstruction efficiency, evalu- 
ated with control samples of r decays. A small correc- 
tion is also included due to a known data/MC difference 
in tracking efficiency. Other sources of systematic un- 
certainty considered include: the number of BB pairs 
in the data sample; the limited size of the signal MC 
samples; the PID efficiency, which is evaluated with con- 
trol samples of A — ► pn~ , D* + — ► L>°(A"~7r + )7r + , and 
4> — > K + K~ decays; possible differences in AE resolution 
between data and MC, which are estimated with control 
samples of B — > DDK decays; charmed baryon branch- 
ing ratios relative to the control modes [l[ ; the A branch- 
ing fraction [l[ ; the presence of intermediate resonances 
in the charmed baryon decay and possible structure in 
the 3-body S-meson decays; and the assumption that 
B(T(4S) -* B°B°) = B(Y(4S) -► B+B~) = 0.5. For fit 
parameters which are fixed to values from fits to the sig- 
nal MC, we vary the value by the uncertainty and take 
the largest change as a systematic uncertainty. Dividing 
out the absolute A + branching fraction also introduces a 
large systematic uncertainty, which we quote separately. 

To investigate whether the three-body mode B~ — > 
A + A~K~ contains intermediate resonances, we exam- 
ine the Dalitz plot structure of candidates in the sig- 
nal region (toes > 5.27GeV/c 2 ), shown in Fig. 03 After 
taking into account the expected background (estimated 
from the toes sidebands), the A + K~ mass spectrum of 
the data is inconsistent with a phase-space distribution 
(% 2 probability of 1.5 x 10~ 7 ). Fitting the data with a 
single, non-relativistic Breit-Wigner lineshape convolved 
with a Gaussian function for experimental resolution, 
we obtain m = 2931 ± 3(stat) ± 5(syst) MeV/c 2 and 
r = 36 ± 7(stat) ± ll(syst) McV. We do not see any such 
structure in the toes sideband region. This description 



TABLE II: Summary of relative systematic uncertainties (%) 
on the branching fractions. The uncertainty on the A + 
branching fraction is 26% and is quoted separately. 



Source 


A+A-K- 


S°A- 


S+A- 


A+A-K 


Tracking efficiency 


9.9 


10.0 


11.4 


11.4 


B counting 


1.1 


1.1 


1.1 


1.1 


MC sample size 


0.8 


1.6 


2.4 


1.5 


PID efficiency 


4.6 


3.5 


3.0 


4.0 


AE resolution 


3.0 


3.0 


3.0 


3.0 


Intermediate BFs 


3.4 


6.9 


0.8 


0.1 


At -> pK~TY + Dalitz 


2.9 


1.8 


1.8 


3.6 


B -> A+A~K Dalitz 


6.9 






4.2 


T{AS) BF 


3.0 


3.0 


3.0 


3.0 


Fit related 


2.0 


1.4 


3.5 


2.5 


Total 


14.5 


13.7 


13.4 


14.3 



is in good agreement with the data (\ 2 probability of 
22%) and could be interpreted as a single resonance 
with those parameters, though a more complicated ex- 
planation (e.g. two narrow resonances in close proxim- 
ity) cannot be excluded. Due to the limited statistics, the 
hclicity angle distribution does not distinguish between 
spin hypotheses. 

In summary, we have studied S-meson decays to 
charmed baryon pairs in four decay modes using a sam- 
ple of 230 million T(4S) — > BB events. The branching 
fraction of B~ — > A + A~K~ is found to be larger than 
the previous measurement [Io| and is comparable to the 
0(1O -3 ) branching fraction predicted for two-body de- 
cays to a pair of charmed baryons. The other results 
are consistent with the previous values 0, E3- The 
data in the Dalitz plot and two-body mass projections 
of B~ — ► A + A~ K~ are inconsistent with a phase-space 
distribution and suggest the presence of a resonance 
in the decay. 
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FIG. 3: Reconstructed B~ —* A+A"^ - candidates in the signal region (?ties > 5.27GeV/c 2 , Ai? < 22 MeV), shown as (a) the 
Dalitz plot, (b) the A^K~ invariant mass distribution, and (c) the A+A~ invariant mass distribution. Data from the signal 
region are shown as black points. Signal events from a phase-space simulation are shown as small grey points in (a) and as a 
histogram in (b) and (c). Data from the sideband region 5.20 < tues < 5.26 GeV/c 2 are shown as a shaded histogram in (b) 
and (c), normalized according to the expected background yield in the signal region. The masses of the B-meson candidates 
are not constrained. 
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